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Abstract 

The first measurement of the elementary process fj,~p — > v^wy is reported. 
A photon pair spectrometer was used to measure the partial branching ratio 
(2.10 ± 0.22) x 10~ 8 for photons of k > 60 MeV. The value of the weak 
pseudoscalar coupling constant determined from the partial branching ratio 
is g p ( q 2 = -0.88m 2 ) = (9.8 ± 0.7 ± 0.3) • g a (0), where the first error is the 
quadrature sum of statistical and systematic uncertainties and the second 



error is due to the uncertainty in \ op , the decay rate of the ortho to para 

pfip molecule. This value of g p is ~1.5 times the prediction of PCAC and 

pion-pole dominance. 
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In semi-leptonic weak interactions the strong force can, in general, induce four couplings 
in addition to the usual vector (g v ) and axial vector (g a ) couplings: weak magnetism (g m ), 
pseudoscalar (g p ), scalar (g s ), and tensor (g t ) ]]]. G-parity invariance predicts null values 
for g s and g t , in agreement with experiment [0-0. The CVC prediction for g m [|J also agrees 
with experiment |5|]. Theoretical calculations of 3% precision are available for g p , but no 
measurement approaching such accuracy has been reported. 

It is difficult to measure g p in beta decay because of the small value of the 4-momentum 
transfer, q. In muon capture the pseudoscalar coupling contribution is larger. Capture on the 
proton allows a reliable calculation since this case is free of nuclear structure uncertainties. 
Assuming partial conservation of the axial current (PCAC) and pion-pole dominance of the 
pseudoscalar coupling, g p can be related to g a (0) 0. In Fearing's notation J8j, 

g p (q 2 = 6.77g a (0)^— 2 (1) 

(7774 — q ) 

where m n is the charged-pion mass and m M is the muon mass. In ordinary muon capture 
(OMC), fi~p — > u^n, q 2 = — 0.88m^, which is far from the pion pole. Thus g a dominates 
g p and a rate measurement of 4% precision || yields a 40% uncertainty in g p . Combining 
many different OMC results yields only 25% precision || for g p . 

In radiative muon capture (RMC), fi~p — > v^wy, q 2 can be much closer to the pion 
pole, since q 2 = at the maximum photon energy (k ~ 100 MeV). The g p contribution 
to the RMC amplitude is then more than 3 times its contribution to the OMC amplitude. 
Offsetting this long known, promising sensitivity to g p is the fact that, in hydrogen, the 
partial branching ratio is ~ 10~ 8 . Due to this small probability and the many potentially 
large background sources, no previous measurement has been attempted. 

The experiment was performed using the M9A muon beam at the TRIUMF cyclotron. 
Negative muons at 63 MeV/c were selected by an rf separator, counted by four beam scin- 
tillators, and stopped in a liquid hydrogen target. The \i~ stop rate was ~6.5xl0 5 s _1 and 
the 7T~ contamination was vr/yU < 2xl0 -4 . Photons emerging from the target were detected 
by 7 ~* e + e~ conversion in a cylindrical lead sheet and analyzed by tracking the e + and 
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e~ paths in cylindrical drift and wire chambers. An axial magnetic field allowed the de- 
termination of the photon energy from the curvature of the e + and e~ paths. Concentric 
segmented cylinders of plastic scintillators just inside the lead converter radius (the A, A', 
and B rings), just outside the lead converter radius (the C ring), and just outside the drift 
chamber outer radius (the D ring) provided the photon trigger via the logic A ■ A' ■ B ■ C • D. 
Lastly, two layers of plastic scintillators and drift chambers surrounding the magnet were 
used to identify cosmic-ray backgrounds. For more details on the experimental setup, see 
Ref. |U|. 

Even a small tt~ contamination in the beam could give rise to a potentially dangerous 
background, since each 7r~ produces photons from the charge exchange or capture reactions: 

n-+p^Ti + n (60.7%), (2) 

->7 + 7 (k ~ 55 - 83 MeV), 
n-+p^-f + n (39.3%, A; ~ 129 MeV). (3) 

Although these pion-induced photons were ~ 10 4 times more plentiful than RMC photons, 
they were prompt and could be vetoed by a timing cut. They were, however, useful for 
measurement of the acceptance and quantitative comparison to Monte Carlo modeling of 
the detector. For this measurement the beam was periodically tuned for 81 MeV/c tt~ , which 
stopped at the same place in the target as the 63 MeV/c \i~ . The pion-induced photons 
were subjected to the same geometry cuts as were the photons from /J,~p, but not to the 
timing cuts. Excellent agreement with Monte Carlo was obtained |[L0| , especially in the 7r° 
region which overlaps the RMC region, providing an energy calibration to ~ 0.25 MeV. 

The spectrometer energy resolution, ~ 11% FWHM for photons, was sufficient to elim- 
inate most of the background for k > 60 MeV from jjT decay with internal or external 
bremsstrahlung, jjT — > e~V~ e v^. These photons have a maximum energy of m M /2, which 
prevents study of RMC at lower energies. Furthermore, fi — > eVwy events could leak from 
the region A; < 60 MeV into the region k > 60 MeV due to the high-energy tail of the 
spectrometer resolution function. To determine this high-energy tail background the photon 
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spectrum from n + stops, where decay is present but capture is absent, was measured. Com- 
parison of the normalized fi + and fi~ photon spectra, accounting for in-flight annihilation 
which contributes to the former but not the latter, yielded a (17.2 ± 2.5)% high-energy tail 
background in the RMC spectrum. 

Some fraction of the cosmic-ray induced events were not identified in the cosmic-ray 
detectors surrounding the spectrometer magnet, and thus survived the software veto. The 
residual background was measured under experimental conditions when the beam was off. 
This and any proton "beam related" background was measured when the secondary beam 
was diverted to other channels. The combined background photon energy spectrum, nor- 
malized by time and integrated proton current, agrees with the spectrum, 100 < k < 200 
MeV, produced by the \x~ beam. This energy range is above the RMC range and allowed a 
reliable extrapolation and subtraction of these backgrounds, (10 ±4)%, in the RMC region. 

Muon capture occurs on the Au flask containing the liquid H and on the surrounding Ag 
heat shields much more rapidly than on H. Au was chosen for the flask material because it 
can be made into a very pure thin-walled container. Also, the disappearance time for fi~ 
in Au is 73 ns, in Ag 89 ns, while in H it is 2195 ns [IT]. Events which occurred within 



365 ns (5t^ u ) following the last /i~ stop were "blanked" in software. A fit to the unblanked 
time spectrum for k > 60 MeV enabled a determination of the part of this background 
which occured more than 365 ns after a stop. The net background from the Au and Ag was 
(10 ±4)%. 

Impurities of deuterium or heavy elements in the liquid hydrogen could have been a 
serious source of backgrounds. In liquid the \i~ is rapidly captured into the singlet atomic 
state. This small neutral system quickly forms a pfip molecule, usually in the ortho state, 
which can decay (A op = (4.1 ± 1.4) x 10 4 s _1 |J) to the para state. However, if there are 
contaminants the yT can transfer from the proton to form a /iA atom. Heavy impurities 
were reduced to <C 10 -9 by thorough bakeout and pumping, and by passing the gas through 
a palladium filter [HJ] . If deuterium is present a fid atom can be formed, followed by a pfid 
molecule. The latter can form a /i 3 He system via muon-induced fusion. Natural H contains 



> 100 ppm deuterium, and since capture on 3 He occurs at a much higher rate than on H, 
the background from 3 He RMC would be larger than the H signal. Our liquid was obtained 
from the electrolysis of two samples of deuterium-depleted H2O. For half of the data-taking 



it contained 1.4 ± 0.2 ppm deuterium ||13|| , while for the other half this contamination was 
< 0.1 ppm. A (1.1 ± 0.3)% 3 He RMC background was determined from a short run with 
natural H. 

Other background sources were explored. These include: 1) tt~ and /i~ capture in the 
Pb beam collimators without hits in the beam counters. As these collimators were located 
quite far from the target such events were vetoed by geometry cuts. 2) Beam tt~ could decay 
in flight to high-energy which would capture on the Au/Cu backplate of the target. The 
beam simulation which reproduced the observed stopping distribution indicated that the 
number of such events was small and a geometry cut further reduced them. 3) /i~ which 
stopped in the veto scintillators would capture on C which has a disappearance time similar 
to H. The stopping simulation indicated this background was < 1%. 4) Before forming a 
pfip molecule the fip atom could diffuse into the Au and capture after a long time, thereby 
defeating the blank cut. The beam simulation showed very few stops within 2 mm of the Au 
flask. The mean free path of such atoms is <C 1 mm and thus this background is negligible. 
5) Multiple single tracks in near-coincidence could imitate high-energy pairs. However, all 
trigger counters were connected to multihit TDCs, which allowed a straightforward veto of 
these spurious events. The net contribution from all these background sources was (3 ±2)%. 

Data taking took place over a 3.5 year period. The total number of fi~ stops examined 
was 3.6 x 10 12 . The uncut photon energy spectrum is shown in Fig. 1. The peaks are 
from bremsstrahlung and pion-induced reactions. The shaded region shows the spectrum 
after application of the prompt cut, indicating elimination of pion events. The number 
of RMC photons observed (k > 60 MeV), after all cuts and background subtractions, is 
N = 279±26. 

The open histogram in Fig. 2 shows the photon spectrum after geometry, timing, cosmic 
ray, and energy (k > 60 MeV) cuts were imposed. The overlayed dark shaded region shows 
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the remaining bremsstrahlung background, obtained from the normalized /i + spectrum (with 
annihilation- in- flight removed). The light shaded region shows the photon spectrum after 
the residual Au, Ag, 3 He, cosmic ray, "beam-related" and bremsstrahlung backgrounds have 
been removed. This represents the first observation of a signal from RMC on hydrogen. 

In order to extract g p from the data, a calculation of the RMC branching ratio _R 7 and 
photon energy spectrum in terms of g p was made using the theory of Fearing and Beder 
This is a relativistic approach based on tree-level Feynman diagrams. Monte Carlo 



events were generated from the theoretical photon spectrum for k>50 MeV, t >365 ns, and 
several values of g p , and then analyzed like experimental events. The q 2 dependence of g p 
assumed was that of Eq. 1. The values of the partial branching ratio and g p were obtained 
when the number of Monte Carlo RMC events iV was equal to the number of experimental 
RMC events JV . 

The number of Monte Carlo events generated was 

N = RyKN^PFa. (4) 

is the number of beam counter coincidences (stops), corrected for multiple /i~ in a beam 
burst. e cr is the efficiency of cuts to eliminate multiple single particles and cosmics (these 
cuts were not applied to pion-induced photons from which the acceptance was obtained, nor 
to simulated photons). It is 0.93 ± 0.01 and was obtained by noting the effect of these cuts 
on fi~ beam events which were rejected by the prompt cut, i.e. 7r~-induced events which 
came in the fi~ beam. Pileup muons blank a fraction 1 — P of real and background events 
which are captured more than 365 ns after a stop, where P = 0.750 ± 0.011 was calculated 
from knowledge of the timing logic and the incident beam rate. This calculation reproduced 
the ratio of events before and after the blank cut, after all other cuts had been imposed. 
F a is the ratio of the measured acceptance in each run period to the simulated acceptance. 
The acceptance was independent of beam rate and within a running period was constant to 
within a few percent, but there were variations of up to 10% between running periods. 

The factor k = 0.88 ± 0.01 is a run-independent product of several small, measured or 



calculated corrections. The largest, 0.926±0.005, is the fraction of muons in that stopped 
in the liquid. It was measured by replacing the liquid with a scintillator of appropriate shape 
and counting coincidences with the beam counters; the measured value agreed with Monte 
Carlo predictions of the beam. Other factors in k, each of which had less than a 2% effect, 
included small measured inefficiencies in the beam and trigger counters, random vetoing in 
the trigger due to the singles rates in the A, A', B veto scintillators, and muon miscounting 
related to electrons in the beam or from the target. 

Due to the spin- dependence of the weak interaction, the extraction of g p from the data 
depends on knowledge of the relative fraction of muons in the ortho and para pp,p states. 
Thus the Monte Carlo, which includes all the muonic atomic and molecular processes, was 
run for various values of X op . The sensitivity to the rates of the other atomic and molecular 
processes was negligible. 

Our measured partial branching ratio for photons of k > 60 MeV is (2.10 ± 0.20 ± 
0.09) x 10~ 8 , where the first error is statistical and the second systematic. This branching 
ratio arises from the mixture of muonic states relevant in the experimental time window 
(t > 365 ns): 6.1% in the singlet [ip state, 85.4% in the ortho p\ip state and 8.5% in the 
para pup state. The shape of the energy spectrum in the experimentally-accessible region 
is in good agreement with the theoretical prediction (see Fig. 2). The results for g p are 
shown in Fig. 3, where g p /g a is plotted versus X op . The expected theoretical result (Eq. 1) 
is 6.77 • g a {0), at q 2 = -0.88m;;. RMC yields g p /g a {0) = 9.8 ± 0.7 ± 0.3, a factor ~1.5 times 
the expected value, if the experimental value X op = (4.1 ± 1.4) x 10 4 s _1 is used. The first 
error includes statistical and systematic errors added in quadrature, while the second is due 
to the error in A op . The present result for g p is independent of energy cut, k > 60 MeV, 
blank time cut, t b > 365 ns, and is not very sensitive to \ op . 

In Fig. 3, "Saclay" refers to the most recent and accurate OMC measurement ||, while 
"World" refers to the world average of OMC (including older bubble chamber) data ||. 
It should be noted that this average was obtained from many different experiments with 
widely varying conditions and results. Within large errors, OMC is in agreement with the 
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expected value of g p if one adopts the measured value of A op . However, the OMC results 
are generally more sensitive to X op than the present measurement, especially the "Saclay" 
result, and the extracted value of g p would decrease significantly if the theoretical value 
\ op = (7.1 ± 1.2) x 10V 1 [0| is correct. 

In conclusion, we have made the first measurement of the elementary radiative muon 
capture process, and determined a partial branching ratio of (2.10±0.22) x 10~ 8 , for photons 
of k > 60 MeV. We have assumed that the measured values of X op and other /zH parameters 
are correct, the q 2 dependence of g p is given by (1), the i? 7 and energy spectrum for RMC 



on H are correctly predicted by a relativistic perturbation theory calculation |14| based on 
tree-level Feynman diagrams, and that a large background has not eluded detection in this 
experiment. Based on these assumptions, a substantial deviation of g p from the predicted 
value is obtained. 

We wish to thank the TRIUMF hydrogen target group, in particular A. Morgan and P. 
Burrill for their excellent work on our target. We are grateful to H. Fearing for access to 
his RMC computer code and helpful discussions. This experiment was supported by the 
NSERC and NRC of Canada, the NSF of the USA, the PSI of Switzerland, and the ARC 
of Australia. 
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FIGURES 
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FIG. 1. Photon energy spectrum before cuts are applied; both ji~ and tt~ induced events are 
observed. Shaded region: prompt cuts applied, pion events are eliminated. 
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FIG. 2. Open histogram: Photon energy spectrum after all cuts were applied. Shaded his- 
togram: final spectrum after background subtraction. Dark shaded overlay: normalized /u + (an- 
nihilation-in-flight removed) bremsstrahlung background spectrum. Curve: theoretical spectrum 
[14] for the best-fit value of g p /g a - 
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FIG. 3. gp(q 2 )/g a (Q)i evaluated at q 2 = — 0.88m^, vs. the ortho to para transition rate A, 
OMC and RMC on hydrogen (see discussion in text). 
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